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Graphene two-dimensional nature combined with today lithography allows to achieve nanoelec- 
tronics devices smaller than the Dirac electrons wavelength. Here we show that in these graphene 
subwavelength nanodevices the electronic quantum transport properties present deep analogies with 
classical phenomena of subwavelength optics. By introducing the concept of electronic diffraction 
barrier to represent the effect of constrictions, we can easily describe the rich transport physics in 
a wealth of nanodevices: from Bethe and Kirchhoff diffraction in graphene slits, to Fabry-Perot 
interference oscillations in nanoribbons. The same concept applies to graphene quantum dots and 
gives new insigth into recent experiments [T] on these systems. 



PACS numbers: 72.10.-d, 72.10.Fk, 42.25.-p, 42.25.Fx 

Analogies play a preminent role in physics. They al- 
low the transfer of notions and concepts from one field to 
another, thus leading to deeper insights and advances in 
both fields. In particular, attempts to explore an anal- 
ogy between the fields of quantum transport and optics 
have been suggested in the past, both in the direction 
of optics, for instance in order to understand coherent 
multiple scattering of light 2 , and also in the reverse 
direction pj]. Optical concepts have also been exploited 
to explain different properties in graphene [4]. In this 
work we develop the analogy between quantum transport 
and optics beyond those attempts and bridge coherent 
electronic transport with sub-wavelength optics. Indeed 
graphene [5j|6j|7], either in the exfoliated or in the epi- 
taxial form, has extraordinary coherence properties. The 
electron mean free-path can be of the order of a micron 
at room temperature and the Dirac electrons wavelength 
can be up to 100 nm or even larger [8]. From a techno- 
logical point of view, graphene can be manipulated by 
ordinary lithography to realize devices as small as few 
nanometers [5] thus achieving a sub-wavelength regime. 
This suggests that an analogy with sub-wavelength op- 
tics can reveal fruitful in order to have a better insight 
into graphene transport properties. We believe that this 
approach could ease the task to find a solution to the 
formidable quantum transport problem of calculating the 
conductance characteristics in nanodevices. 

We will show that, by introducing the hybrid 
"electronics-optics" concept of electronic diffraction bar- 
rier to represent constrictions at contacts, we can de- 
scribe in an easy way the rich transport physics in a 
wealth of sub- wavelength graphene nanodevices: from 
the simplest systems such as slits between two half 
graphene planes, to graphene nanoribbons or quantum 
dots sandwiched in between semi-infinite graphene sheets 
as leads [10] [HI El EGS O EH [16] . Our methodology 
relies on an exact numerical calculation of the conduc- 
tance within a tight-binding |17j and a Landauer quan- 



tum transport formalisms [3]. The novelty here intro- 
duced is that the Landauer equation is solved via a new 
recursive numerical algorithm [18 . 

Our results give new insights on the conductance char- 
acteristics of experimentally synthctized graphene quan- 
tum dots [H [19] where a chaotic Dirac billiard behaviour 
has been recently observed pQ. The diffraction barrier 
concept explains electron localization and allows us to 
provide a new interpretation of experimentally observed 
features observed into this system as due to diffraction 
phenomena occurring at the divergent wavelength of elec- 
trons at the Dirac point. We believe that the analogy 
can be further pushed to play a fruitful and key role 
in the design of new graphene nanodevices exhibiting 
other original and interesting characteristics. We also 
point out that the present findings are not only specific 
of graphene Dirac electrons but could be extended to 
other sub-wavelength devices made by 2D electron gas 
such as GaAs/AlAs or in molecular electronics. 

Conductance of graphene slits and diffraction by an 
aperture - The archetype experiment showing an optical 
diffraction behaviour is the transmission through a slit 
or hole. We considered slits made into a graphene sheet 
(Fig. [T]) [24 and calculated the conductance. We took 
into account several slits differing by the width VP, from 
the thinnest consisting only of a single motif i.e. a single 
hexagon. We introduce the wavelength A of the incident 
Dirac electrons, related to the channel energy E via A = 
2W _ kv^ where j s Dirac electrons wavevector, E = 
hkvp is the energy, linearly dispersed in graphene with a 
Fermi velocity of vp ~ 10 6 ms" 1 [B]. 

In Fig. [T] we report the conductance of the slits as a 
function of VP/A. We first remark (Fig. [I] inset) an uni- 
versal scaling behaviour independent of the slit width 
W and depending only on the ratio VP/A. This scaling 
breaks down when the wavelength is A < 4 nm, that is 
the energy is E > 1 eV. Indeed, as shown below, this 
scaling is intimately related to the scaling properties of 
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FIG. 1: Diffraction of electrons through a graphene slit. Left: 
schematic view of the geometry of a graphene slit. Right: 
universal scaling law of the graphene slit conductance as a 
function of W/X. This shows the Bethe-like X^$> W quadratic 
regime and the Kirchhoff-like A <C W linear regime (see text). 
The inset shows the conductance for slits of widths W — 0.8p 
nm with p — 1, 2, 3... 



the Dirac equation which in graphene is valid only in the 
low energy or long wavelength limit. 

We then notice that in the limit W/X <C 1, the con- 
ductance is quadratic, with the conductance g(W/X) ~ 
50(e 2 /h)(W/X) 2 , while in the opposite limit W/X > 1, 
the conductance turns out to be linear with g(W/X) ~ 
5(e 2 /h)(W/X). The crossover between the two regimes 
occurs around W/X ~ 0.2 ~ 0.5. This promptly re- 
minds of an analogy with classical optics and hence of- 
fers an immediate interpretation: the conductance re- 
sponse of a graphene slit is a clear manifestation of a 
diffraction phenomenon. Indeed we can identify two dif- 
ferent diffraction regimes: for wavelengths X^$> W, much 
larger than the aperture, we observe a Bethe-like diffrac- 
tion regime [10] with a slit transmitted energy propor- 
tional to the square of the aperture. On the other hand, 
for A <C W one has a Kirchhoff-like diffraction regime, 
the transmitted energy being proportional to the aper- 
ture. Note that weak-localization effects in quantum bil- 
lards can also reveal diffraction effects in the Kirchoff 
regime at the lead mouth [5D]. In the short wavelength 
regime the phenomenon of diffraction is a perturbation 
and a semi-classical description of the electrons as ballis- 
tic wavepackets is applicable, thus leading to a transmit- 
ted current proportional to the section of the slit W. By 
contrast, in the subwavelength regime the transmission 
is lower than given by the semi-classical picture and can 
be viewed as a tuneling process. 

In the ordinary optical diffraction of the slit transmit- 
ted intensity, an universal scaling behaviour is observed 
due to the scaling invariance of Maxwell equations. Here 
it stems from the fundamental scaling invariance of the 
graphene 2D Dirac equation. Indeed if a spinor ip(r) sat- 
isfies the Dirac equation for the energy E and wavelength 
A then the spinor ip(r) = ip(xr) satisfies the Dirac equa- 
tion for the energy E = E/x and wavelength xX. This 
leads to the fact that the conductance undergoes a scale 



invariance g(xX,xW) = g(X, W). As a consequence the 
conductance only depends on a reduced argument W / A, 
g(X,W) = g(W/X). [25] , Let us note that a scaling rela- 
tion can also exist for other geometries. In that case the 
conductance of the circuit will depend only on the ratio 
between the characteristic lengthes of the circuit and the 
wavelength A. 

Nanoribbons as Fabry-Perot interferometers and sub- 
wavelength waveguides - Here we study the quantum 
transport response of finite-length metallic graphene 
nanoribbons and show that they present an oscillating 
response, typical of optical interferometers, and behave 
as subwavelength waveguides. We consider two ribbons, 
a zigzag and an armchair both chosen with a metallic 
character and hence an available conductance channel 
at the Dirac energy [25]. The geometry is presented in 
Fig. [2] and we have taken a length of 6 nm for the arm- 
chair nanoribbon and 3 nm for the zigzag (see Fig. [3] 
for the detailed geometry of each junction). The calcu- 
lated conductance for both ribbons is shown in Fig. [2] 
where we again observe a conductance going to at 
the Dirac point. However in this case, the most evident 
features are large amplitude oscillations in the conduc- 
tance, from maxima values of 1 to minima disposed along 
an envelope function. These systems behave as Fabry- 
Perot cavities, exactly as in subwavelength optical metal- 
lic waveguides [TH [T31 ES] • To demonstrate this we show 
that the nanoribbons calculated oscillating conductance 
distributes like the Fabry-Perot standard transmittancc 
which is the Airy function: 



Tfp(E) = 



1 + F(E) suS(4>(E)/2) 



(1) 



where cj>{E) = 2k(E)L + 24>{E) is the phase difference 
after one loop in the nanoribbon. k(E) is the wavevector 
of the Bloch state in the infinite ribbon with energy E. 
4>(E) is the phase factor acquired at each reflection. L 
is the length of the Fabry-Perot interferometer. F(E) = 
AR(E)/(l - R(E)) 2 is the finesse coefficient and R(E) = 
1 — T(E) is the reflection coefficient at each end of the 
nanoribbon, that is at the junction. 

The Airy function presents maxima T^ LX — 1 equal to 
1 when cj)/2 is an integer m multiple of it. For sufficiently 
large L the phase 4>{E) = 2k(E)L + 2(f>(E) varies rapidly 
with energy, as compared to F(E); and the minima occur 
essentially when sin2(</>/2) is maximum, that is at <j)/2 = 
ran + tt/2, and envelope along the function 



FP \&) 



1 ^ (l-R(E)Y 
l + F(E) (l + R(E)) 2 



(2) 



which is independent of the length of the ribbon. We 
indeed checked that the envelope functions of the con- 
ductance minima, as expected from Eq. pi), is indepen- 
dent of the ribbon length. In particular the minimum of 
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FIG. 2: Graphene nanoribbons as Fabry-Perot interferom- 
eters. Left: schematic view of nanoribbons contacted to 
graphene half-planes. Right: Fabry-Perot oscillations inf 
the conductance of armchair (upper panel) and zigzag (lower 
panel) nanoribbons. The dashed line which almost perfectly 
envelops along the oscillations minima, is the plot of Eq. Q 
with R(E) taken from Fig. [3] 

the transmission tends to zero at the Dirac energy, which 
indicates that R(E) tends to 1 at zero energy. 

Another feature of the Fabry-Perot resonnances con- 
cerns their width. The peaks full width at half maximum 
is ^fwhm = H% = where 5<j) = 2(1 - R)/VR 

and v(E) is the group velocity at energy E. The velocity 
is bounded, such as the peaks look thin when R is close 
to 1, and broaden when R — > 0. In our case when E — > 
the peaks seem to have negligible width, indicating again 
that R — > 1, and broaden for the higher energies, indi- 
cating lower values of R. 

The relative variation of the conductance (i.e. the ratio 
between the maximum and minimum conductance be- 
tween two peaks) is also maximum close to the Dirac 
energy, due to the nearly perfect confinement at this en- 
ergy (R(E) ~ 1). We consider in the following that reso- 
nances are well defined if the maximum to minimum ratio 
is greater than 2. In our calculations this corresponds to 
a criterion X/W > 3 ~ 5 where W is the width of the 
nanoribbon. 

The concept of electronic diffraction barrier at a con- 
striction - Fabry-Perot oscillations occur due to the re- 
flection at the junctions/ends of the nanoribbons. We 
now analyse this reflection and the response of the ele- 
mentary junction. We consider the same metallic zigzag 
and armchair ribbons as before, but take them as semi- 
infinite (Fig. [3]) . If the system were ballistic (infinite rib- 
bon, no contact resistance and no reflexion), we would 
expect to observe a conductance equal to 1 quantum of 
conductance (2e2/h) (dot-dashed line in Fig. |3| in the 
range of the Dirac energy. Instead, the calculated trans- 
mittance T(E) through the junction (continous line in 
Fig. [3]) is observed to drop to zero at the Dirac energy, no 
matter the chirality, exactly like in slits and the nanorib- 
bons. The only effect of the ribbon electronic structure is 
to produce a little difference in the characteristics, with 
the zigzag conductance showing a more cusp-like feature 
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FIG. 3: Electronic diffraction barrier at a contact. Left: the 
contact geometries of the armchair (upper panel) and zigzag 
(lower panel) semi-infinite ribbons coupled to graphene half- 
planes. Right: Conductance g of the device represented on the 
corresponding left side, g = T2e 2 /h with T the transmittance 
of the diffraction barrier. 



at the Dirac energy. 

Since the associated resistance is independent of the 
ribbon length, this behaviour is a genuine manifesta- 
tion of the contact resistance at the junction [26 . The 
junction conductance we have here calculated represents 
the characteristics of the elementary electronic diffrac- 
tion barrier which is at the basis of the universal be- 
haviour observed in all the devices hereby studied, slits, 
nanoribbons or quantum dots. Indeed coming back to the 
nanoribbon case, if in Fig. [2] we plot (1 — R(E)) 2 /(l + 
R(E)) 2 (dashed curve) with R(E) = 1 - T{E) as ex- 
tracted from Fig. [3] for the elementary junction, we find 
that the curve turns out to be the envelope of the Fabry- 
Perot oscillations minima in the nanoribbon conduc- 
tance, as expected from Eq. ([2]). This exemplifies the 
concept of electronic diffraction barrier at a constriction 
and gives a quantitative example of such a barrier. 

Conductance of a quantum dot - We now turn to more 
complicated structures and consider a standard nano- 
electronics quantum dot (Fig. [4} consisting of a pur- 
posely irregular shape graphene nanostructure contacted 
via small apertures to the two half-planes graphene leads. 
Applying again the concept of electronic diffraction bar- 
rier, we expect that the quantum dot is weakly coupled to 
the leads when the energy of the electrons in the graphene 
sheet is close to the Dirac energy, i.e. when their wave- 
length is sufficiently large. We note that the conductance 
of this system (Fig. [4]) presents a maximum around 3 
eV which is the energy of the Van-Hove singularity of 
bulk graphene where the density of states diverges. We 
also note many irregular peaks superimposed to the en- 
velope, clearly showing characteristic resonances of the 
quantum dot to be associated to the particular (irreg- 
ular) shape. The conductance decreases when going to 
the Dirac point, which is a signature of the lowest num- 
ber of available states and of a diffraction barrier at the 
contact constrictions. The maximum to minimum trans- 
mission ratio between adjacent peaks is much larger in 
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FIG. 4: Conductance of the quantum dot, showing well de- 
fined resonances close to the Dirac energy. Inset: geometry 
of the irregular shaped quantum dot coupled to the graphene 
half-planes. 



the vicinity of the Dirac energy than at other energies. 
This is also observed in the nanoribbons studied previ- 
ously. Again this can be explained by the presence of two 
diffraction barriers at the contacts which are nearly com- 
pletely reflecting (R — 1) at the Dirac energy. This leads 
to well defined states within the quantum dot, and thus 
well defined conduction channels through these states. If 
we take as criterion to identify a range with well defined 
resonances the same we already introduced (maximum 
to minimum ratio greater than 2), then we find that the 
critical wavelength is of order of A ~ 4nm for a width 
W = 1.2nm. This is the same range of values of X/W as 
in nanoribbons, i.e. X/W > 3 ~ 5. 

Experimental implications - Our quantum dot com- 
puted characteristics (Fig.[4j looks very similar to the ex- 
perimental result of Ponomarcnko et al. (Ref. [1], Fig. 1). 
In those experiments the relation between voltages V 
(here V is the voltage variation with respect to the Dirac 
point voltage expressed in Volts) and electronic wave- 
length A is A ~ 130nm/y/V [23]. From that one deduces 
that in their experiment the resonances are well defined 
for a wavelength A to aperture W ratio X/W > 3 ~ 5 
in reasonable agreement with the numerical results pre- 
sented here. They reported also the low bias conduc- 
tance of constrictions of 20 nm (Ref. pQ, Fig. S2), similar 
to the slits studied here. First of all, after our study of 
slits the Kirchoff-like (linear) regime is reached at volt- 
ages V > 10 Volts and for higher voltages we predict for 
this 20nm slit a conductance g(W/X) ~ 5(e 2 /h)(W/ X) ~ 
0. 75 (e 2 /h)y/V. This is consistent with the experimental 
results and for instance at 100 Volts from the voltage cor- 
responding to the Dirac point they measure G ~ 5.5e 2 //i 
(Ref. pQ, Fig. S2), while we predict G ~ 7.5e 2 /h, in good 
agreement with the experiment. On the other hand, the 
comparison of our results with their smaller 10 nm con- 
striction is not such good. Other mechanisms should 
influence the electronic transport in this case. This is 
in agreement with the evidence of thermally activated 
transport for the smaller experimental constriction. 

To conclude, our work establishes an important link 



between nanoelectronics and subwavelength optics which 
could be tested in today's technology graphene nanode- 
vices, thank to its remarkable properties. The concept 
of electronic diffraction barrier is central to the under- 
standing of transport properties in the subwavelength 
regime. It allowed us to get insights into the proper- 
ties of graphene junctions, slits, nanoribbons, and quan- 
tum dots and explain their conductance characteristics 
in terms of diffraction and interference phenomena. This 
analogy can be further pushed forward in the design of 
devices with new properties. 
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